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I . Abstract. - We examine the crossover from classical to non-classical critical behaviour in two- 

' ^ ' dimensional systems with a one-component order parameter. Since the degree of universality of 

C ' the corresponding crossover functions is still subject to debate, we try to induce non-universal 

O , effects by adding interactions with a second length scale. Although the crossover functions 

, ^, . clearly depend on the range of the interactions, they turn out to be remarkably robust against 

further variation of the interaction profile. In particular, we find that the earlier observed 

non-monotonic crossover of the effective susceptibility exponent occurs for several qualitatively 

►^ , different shapes of this profile. 
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^\ ' Introduction. - In recent years, there has been a revived interest in the nature of the 

0^ ■ crossover from classical to non-classical (asymptotic) critical behaviour upon approach of the 

critical point. This crossover between two universality classes can be observed in a great variety 

of many-body systems, including pure fluids, polymer mixtures and micellar solutions, and is 

driven by the ratio between the reduced temperature t = {T — Tc)/Tc (where T^ is the critical 

i-j^ . temperature) and a system-dependent parameter G, the Ginzburg number. The dependence 

r^ ' of observables on this ratio is described by so-called crossover functions. Compared to our 

O ■ knowledge of critical exponents, for which very accurate, consistent estimates are available 

J-^ ' from renormalization-group (RG) calculations, series expansions, experiments and numerical 

calculations, the situation for non- asymptotic critical phenomena is not so clear-cut. Most 

theoretical predictions for crossover functions are obtained by means of RG-based methods. 

Examples include the work of NicoU and Bhattacharjee [|l|, who used an RG-matching method 

C^ ' to calculate crossover functions for the one-phase region (T > Tc) to second order in e = 4 — d 

{d denotes the spatial dimensionality), and of Bagnuls and Bervillier pj, who applied massive 

field theory in d = 3. The latter work was then extended to the two-phase region as 

well PI, although only for temperatures relatively close to T,,. Very recently, the approach 

of Ref. I^was used to calculate the full crossover function for the susceptibility exponent 

below Tc |H] . A more phenomenological approach has been taken by Belyakov and Kiselev H , 

who presented a generalization of flrst-order e-expansions. Although the different methods vary 

in mathematical rigour, they all suggest that the crossover functions are universal functions 

of the ratio t/G, under the additional restriction i ^ 0, G — > 0. This limit, which is referred 
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to as critical crossover, implies that one must consider the limit in which the coefficient u of 
the quartic term in the Landau-Ginzburg- Wilson (LGW) Hamiltonian goes to zero (as can 
be realized, e.g., in systems with a diverging interaction range H). Experimental systems 
evidently do not obey these restrictions: Here G is a fixed parameter and the crossover 
functions are obtained by varying i, where it is generally assumed that within the critical 
region, i.e. for t sufficiently small, one still observes a universal crossover. Recent work [[7|-y 
has shown that this assumption is only partially correct. In this paper, we therefore examine 
the role of some of the parameters that might be held responsible for deviations from the 
field-theoretic crossover curves and show what degree of universality one may still expect. 

Anisimov et al. [|0| have suggested that, apart from the correlation length f , an additional 
(mesoscopic) length scale may determine the nature of the crossover behaviour for complex 
fluids. Based upon earlier work [|ll| , a corresponding parametric crossover function was 
proposed and in Refs. [[7|, ^ it was shown that this function of two variables can indeed 
describe the crossover of the susceptibility exponent for several experimental systems displaying 
qualitatively different behaviour. On the other hand, there are quite a number of experimental 
results that can be described in terms of the above-mentioned sm^/e-parameter functions, 
although it should be noted that only few experiments have yielded accurate results for 
the effective exponents, which are defined as the logarithmic derivatives of the crossover 
functions M]. Furthermore, most experiments only partially cover the crossover region, 
which occupies several decades in the reduced temperature. Thanks to recent algorithmic 
developments, numerical methods can circumvent both of these limitations in an efficient way, 
which has already led to several notable results |^,|l^,|l^. In particular, it was demonstrated 
in Ref. H that the crossover function for the effective susceptibility exponent ^'^g (pertaining 
to T > Tc) for three-dimensional systems with a finite interaction range R is steeper than the 
functions presented in Refs. |0,H]. The reason for this discrepancy lies in u not being small 
for small R. Whereas the scale of the crossover is determined by the Ginzburg number, the 
shape of the crossover functions is determined by m [y. This makes it difficult to obtain an 
explicit expression for the crossover. One possibility is to invoke the description of Ref. . A 
(still somewhat phenomenological) fit to this description is indeed possible ||], but a further 
demonstration that the crossover functions depend on more than one parameter is clearly 
desirable. The numerical results presented in Refs. |8||l^,|l^ have been obtained for a block- 
shaped interaction profile (the so-called equivalent-neighbour model), where the interaction 
strength is kept constant within a radius Rm and zero beyond that. As long as the interaction 
has a finite range, different interaction profiles will lead to the same universal properties |^, 
but not necessarily to the same crossover functions. Thus, there is a twofold objective in 
studying the effect of a modified interaction profile. In the first place, we want to study the 
effect of introducing an additional length scale in the block-shaped interaction profile, in order 
to study its influence on the crossover functions. Secondly, this modiflcation of the interaction 
profile allows us to examine, on a more general level, the dependence of crossover functions 
on the shape of the interaction profile and thus to shed some light on the universal nature of 
these functions. 

Simulational aspects and determination of critical properties. - In order to maximize the 
numerical sensitivity to variations in the crossover function, we have restricted ourselves to 
the two-dimensional (2D) case only. An additional attractive aspect of this system is, that the 
effective susceptibility exponent 7^ exhibits a remarkable non-monotonicity p4| . It is well 
possible that such a peculiar property is particularly sensitive to variations in the interaction 
profile. Thus, we have carried out Monte Carlo simulations for a 2D Ising model, defined on 
a square lattice of size L x L with periodic boundary conditions. Each spin interacts with 
a strength Ki with all its neighbours within a distance r < Ri (domain Di) and with a 
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Fig. 1. - The interaction profile studied in this work. Each spin interacts with a ferromagnetic 
coupling Ki with all its neighbours within a distance Ri and with a coupling K2 with those at a 
distance between Ri and i?2. 



strength K2 with aU its neighbours within a distance Ri < r < R2 (domain 02). This means 
that the block-shaped interaction profile of Ref. M has been generalized to the double-blocked 
case depicted in Fig. n. The strength ratio K1/K2 is denoted by the parameter a, which 
throughout this work is supposed to be greater than unity. In order to suppress lattice effects, 
all range dependences are expressed in terms of the effective interaction range R, defined by 
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which for our interaction profile reduces to i?^ = (a ^ 
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azi + Z2, with Zi the number of neighbours in domain Di. The strength ratio was chosen as 
a — 16, in order to create a strong asymmetry between the two domains. The value of R2 



Table I. - Some properties of the interaction profiles studied in this work. The systems are listed in 
order of increasing R, except for the last one, which is the only profile for which Ri < -Rmin (Ri ~ 27 
is the borderline case corresponding to the minimum in R). 
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0.00872817 (2) 
0.002918960 (7) 
0.000836683 (2) 
0.000598920 (2) 
0.0003934150 (11) 
0.0002775163 (9) 
0.0002140278 (7) 
0.0001777575 (4) 
0.001683755 (5) 



0.8559 (3) 
0.8566 (8) 
0.8551 (8) 
0.856 (2) 
0.858 (3) 
0.856 (5) 
0.856 (5) 
0.850 (4) 
0.859 (4) 



1.8744 (9) 
1.8739 (14) 
1.8738 (13) 
1.873 (2) 

1.872 (6) 

1.873 (3) 
1.878 (4) 
1.867 (5) 
1.873 (5) 
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was kept fixed at yliO. In the finite-size scaling analyses, the minimum system size has to be 
of the order of i?| and a maximum linear system size L — 1000 thus implies that the results 
cover a factor 7 in L. The effective interaction range R was then varied by varying Ri. Both 
for i?i — !■ and for i?i — > i?2, R will take its maximum value (which in the continuum limit 
approaches i?2/V2) and it will reach a minimum at i?i — i?inin- In the continuum limit, the 
corresponding effective range is R^ — R^^^ = R\/{^/a + 1). Although the same values for R 
can be reached with i?i < i^min and i?i > i?min, it should be noted that the two cases greatly 
differ. For example, for i?| — 140 and a = 16 (where iJ^jj^ — 28.06 is reached for zi = 88, i.e. 
26 < Rf < 28) one may obtain i?^ ss 50 by choosing either i?^ = 4 or i?^ = 93, but in the 
former case Di contains 12 out of 436 interacting neighbours, compared to 292 out of 436 in the 
latter case. This means that the integrated coupling ratio azi/z2, which indicates the relative 
contribution of the two domains to the total integrated coupling, is 0.45 in the first case and 
32 in the second case. So, in combination with the original block-shaped profile, we can realize 
three qualitatively different interaction profiles and study the dependence of the crossover curve 
on the profile. Table m lists some properties of the interaction profiles considered in this work. 
The three profiles with i?2 < 140 were added in order to reach very small effective interaction 
ranges as well. For each choice, we have carried out extensive simulations using a dedicated 
cluster algorithm for long-range interactions pq| . The critical properties of each individual 
system were determined via finite-size scaling analyses, along the lines described in Ref. m. 
The critical coupling, for which an accurate value is required to attain the proper crossover 
curve, has been determined from the amplitude ratio Q = {m?)'^ / {ra^) and we have obtained 
the magnetic exponent j/h from the absolute magnetization density (|m|) (see Table |). One 



notes that for all systems Q is in good agreement with the 2D Ising value Qi ~ 0.856216 |17 
and yh lies very close to 15/8. This confirms the expectation that all systems belong to the 2D 
Ising universality class. In addition, we have determined the magnetic susceptibility for t < 
from the fluctuation relation x = L'^{{m'^) — (|mp))/fcBT and the critical finite-size amplitudes 



Range dependence of critical properties and analysis of the crossover functions. ~ Figure y 
shows the critical temperatures as a function of the interaction range. All temperatures are 
expressed in units of the critical temperature of the mean-field model, Tc — l/{zcfiKc). 
Since fluctuations are less suppressed when the interaction range decreases, one observes 
that Tc is gradually depressed for smaller R. More importantly, the figure illustrates that 
a definitely non-universal quantity like Tc does not depend on the effective interaction range 
alone: Although for several systems the critical temperature lies on the curve describing Tc{R) 
for the block-shaped interaction profile, the systems with Ri < i?min (-Rf = 4, 27) exhibit 
a clear deviation from this curve. Apparently, the latter systems show (for a given R) the 
greatest deviation from the mean-field model, in the sense that the suppression of fluctuations 
is least efficient. 

In contrast, no such deviations are observed for the critical finite-size amplitude of, e.g.^ the 
absolute magnetization density. This quantity, defined as do = limL^oo i'*^^''(|?Ti|), has an 
asymptotic range dependence proportional to /j(3'^-4yh)/(4-d) g^ j-j- ^^.^g q^^ ^^j-^g^^^ ^j^g systems 
investigated in this work do not only follow this asymptotic law, but also for smaller R agree 
very well with the range dependence found for the block-shaped profile, see Fig. g. This even 
holds for highly asymmetric profiles such as Rf — A]R\ — 140. 

The quantity of central interest, however, is the magnetic susceptibility x below Tc. In Fig. H 
the crossover function for the block-shaped profile is shown as a reference curve. The parameter 
along the horizontal axis is proportional to t/G (G ex _R^2d/(4-d)-j ^^^ ^j^g susceptibility is 
divided by a factor R"^ to obtain a data collapse for different ranges. One clearly observes 
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Fig. 2. - The critical temperatures of the various models as a function of the interaction range. The 
open squares and the curve refer to the block-shaped profile of Ref. H and the black squares to the 
profiles studied in the present paper. The numbers indicate the value of the parameter Rf (with 
R2 = 140 fixed); the indication 25/32 refers to the system with Ri — 25 and R2 = 32. 
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Fig. 3. - Critical finite-size amplitude of the absolute magnetization density. The dashed curve 
indicates the RG expression fitted to the open squares; clearly, it also describes the black squares 
(referring to the interaction profiles of the present work) very well. The numbers are explained in the 
caption of Fig. Q. 



how the solid curve interpolates smoothly, but with a non-monotonic derivative, between the 
Ising asymptote and the classical asymptote. Within the same figure, we have also plotted 
the finite-size data for four different double-block profiles. Several remarks are in order here. 
First, all data have been divided by a range-dependent factor describing the deviation of the 
connected susceptibility from its asymptotic range dependence. This is similar to the difference 
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Fig. 4. - The crossover function for the connected susceptibility below Tc. X = x/C{R), where 
C{R) is a range-dependent correction factor that accounts for the fact that the critical amplitude 
for small R deviates from the asymptotic range dependence. This only introduces a shift along the 
vertical axis. The solid curve indicates the crossover function for the block-shaped interaction profile 
and the dashed lines mark the mean-field ( "MF" ) and the Ising asymptote. The numbers in the key 
refer to the values for R^ and i?2 for each interaction profile. For clarity not all systems listed in 
Table | have been included in the graph, but the omitted data points are fully compatible with those 
shown. 



between the dashed curve and the solid line in Fig. p| as this curve turns out to depend solely on 
the value of R and not on the shape of the interaction profile, it is permissible to use the same 
expression for all systems. For large interaction ranges, the correction factor approaches unity. 
Secondly, at the right-hand side of the graph the data points start to deviate from the reference 
curve. This is caused by the fact that, sufficiently close to Tc, the diverging correlation length 
is truncated by the finite system size. For the systems with R\ — A and R\ = 93 (both with 
i?2 = 140) this happens in the figure at different temperatures, despite the fact that they have 
very similar values for the effective range R. The reason for this is that the data points pertain 
to different system sizes, viz. L = 1000 and L — 300, respectively. The inset shows that for 
the same system size {L — 300) the data points for both systems virtually coincide, even in 
the finite-size regime! Finally, the left-most data points have been corrected for saturation 
effects, which are fully described by mean-field theory, cf. Ref. [|5|. This is merely an optical 
issue: Also the saturated curves (which display a strong decrease of the susceptibility) show no 
dependence on the shape of the interaction profile. The primary message, however, of Fig. 
is that for all interaction profiles the data in the thermodynamic limit perfectly coincide with 
the reference curve for the block-shaped potential. We view this as a strong indication that 
crossover functions possess a considerable degree of universality and conclude that a second 
length scale (-Ri) of the form introduced in this work is insufficient to induce modifications of 
the crossover functions, contrary to some expectations. 



Conclusions. ~ In summary, we have examined the critical properties of two-dimensional 
Ising-like models with an extended interaction range, for several different shapes of the in- 
teraction profile. In addition, we have calculated the crossover function for the susceptibility 
in the low-temperature regime, describing the crossover from classical to non-classical critical 
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behaviour upon approach of the critical point. Although recent work has suggested that this 
function cannot be described in terms of a single parameter, namely the reduced temperature 
divided by the Ginzburg number, we find that it is independent of the precise shape of 
the interaction profile. Irrespective of the presence of an additional length scale or a high 
asymmetry in the interaction profile, all examined systems can be classified according to 
a single additional parameter describing the effective interaction range. In particular, the 
non-monotonic crossover of the effective susceptibility exponent, as found in Ref. (l4|, is 
not a peculiarity of the block-shaped interaction profile, but can be observed for all systems 
studied in the present work, provided that the effective interaction range is sufficiently large. 
Furthermore, a corollary of the results presented here is that the coefficient u in the LGW 
Hamiltonian also appears to have only a weak dependence on the shape of the interaction 
profile. Of course, this still leaves the possibility that other parameters, that still need to be 
identified, have a more pronounced influence on u and thus on the crossover functions. 
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